Decoherence in nonclassical motional states of a trapped ion 
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The decoherence of nonclassical motional states of a trapped 9 Be + ion in a recent experiment is 
investigated theoretically. Sources of decoherence considered here destroy the characteristic coherent 
quantum dynamics of the system but do not cause energy dissipation. Here they are first introduced 
phenomenologically and then described using a microscopic Hamiltonian formulation. Theoretical 
predictions are compared to experimental results. 
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I. INTRODUCTION 

The experiments of Meekhof et al Jj[] have re- 
vealed quantum dynamics characteristic of the Jaynes- 
Cummings type (especially collapses and revivals of exci- 
tation probabilities) Q for the first time in a trapped ion 
system. Stimulated Raman transitions coupled the inter- 
nal states of a trapped 9 Be + ion to its motional states, 
within the Lamb-Dicke limit of tight ion motion confine- 
ment in the trapping potential. The Jaynes-Cummings 
spin-boson Hamiltonian then derives from the coupling of 
the internal electronic states of the ion to the vibrational 
quantum states of motion. 

The characteristic quantum dynamics (collapse and re- 
vival) of the Jaynes-Cummings type interaction for the 
ion motion HH (in the experiment of Meekhof et al [Q , 
an "anti Jaynes-Cummings interaction" for driving the 
first blue sideband) were observed in the population of 
the lower atomic state (P\), which was modelled by the 
phenomcnological form fitting the observation as 



/Jpn cos (2gty/n + l) e 



7n* 



(1) 



Here p n is the initial probability distribution for the mo- 
tional states in the Fock state basis, g is a coupling 
constant between the motional states and atomic states 
(Rabi frequency), and 7„ is a phenomenological damp- 
ing rate. The observed damping rate can be written as 
7n = 70(71 + l) w with v w 0.7 observed in the experi- 
ments of Ref pi] . The damping rate of the nth compo- 
nent is independent of that of different components, so 
that equation (^) implies decoherence without there be- 
ing transitions between the states of different quantum 
numbers (energy relaxation). The conventional sources 
of decoherence, such as spontaneous emission between in- 
ternal atomic states, and population decay of motional 
states, cause transitions between the states of different 
quantum numbers and do not give the decay rate in a 
form which can be written as 7„. There have been sug- 
gestions [U as to the origin of this decoherence with the 
unusual observed value of v, in terms of decoherence of 
the ion motion, decoherence of the ion internal levels, 



and decoherence caused by non-ideal applied fields but 
the situation has not yet been satisfactorily resolved. 

In this paper, we introduce phenomenologically new 
sources of decoherence, which destroy the characteris- 
tic Jaynes-Cummings type dynamics without energy re- 
laxation, by coupling the spin-boson system to a quan- 
tum reservoir The reservoir consists of many- mode 
bosons described by a canonical distribution at tempera- 
ture T and introduces noise to the system. We treat de- 
coherence microscopically using a master equation. The 
master equation coincides with that for stochastic white 
noise in the high temperature limit of the reservoir under 
certain approximations (Markovian approximation and 
ohmic density of states of the reservoir §). The ad- 
vantage of using a quantum boson reservoir is that it 
not only describes phenomenological quantum noise, but 
also gives more microscopic information on the source 
of decoherence, e.g. the noise frequency being responsi- 
ble for decoherence even in the high temperature limit. 
Using this combined approach from two directions (phe- 
nomenological and microscopic) we discuss the origins of 
decoherence in this system. 



II. DESCRIPTION OF THE SYSTEM WITHOUT 
DECOHERENCE 

Before investigating decoherence, we consider the sys- 
tem without decoherence, reviewing how the stimulated 
Raman transitions describe the "anti Jaynes-Cummings" 
interaction (^-||] when the first blue sideband is driven, 
and introducing the dressed states description of the anti 
Jaynes-Cummings system. We note that the first red 
sideband driven case (the Jaynes-Cummings interaction 
case) can be treated just in the same manner, where we 
exchange the two relevant internal atomic levels ||) and 
|t) of the following formulation. 

We consider a system with three internal levels \j) 
(j = 0,1, t) and their motional states \n) (n = 0, 1, ....). 
They are represented by the following Hamiltonian: 



Ha — H ato 



H, 



rib 



(2) 



where 
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Hatom = -fiW0l|l)(l| 

H vib = hw x tfb 



fiW02|T>(T|, 



(3) 
(4) 



with the transition frequency u>oi (^02) between states 
ID (|D) and |0), the creation (annihilation) operator of 
the motional states w (b), and the frequency of the mo- 
tional states uj x . We employ two driving laser beams 
with detuning A, momentum ki (k 2 ) and frequency u>i 
(1V2) which cause dipole transitions between the level \l) 
(ID) and |0). (See Fig. [l]) These beams can be treated 
classically, so the interaction Hamiltonian is 



Hint — —Hoi ■ R- c 



E m e' (krr ^ lf) 



— /Lt 02 • Ke Jio2e 



(5) 



where fi 01 (fi 02 ) is the dipole matrix element between \ l) 
(IT)) and |0). 



lj=0> 




|j=f , n =0 > 



lj=J, , n=0 > 



FIG. 1. Energy levels of the internal states and the mo- 
tional states. 

We apply the rotating wave approximation to the inter- 
action Hamiltonian (||) , transform to the interaction pic- 
ture {Hj nt = e lHot/h H int e- lHot/h ) and expand in terms 
of the motional state quantum numbers. When the blue 
sideband is driven, we have 



= -hJ2goie l{{n ~ m)uJ * +A}t (n|e lkl - r |m) \n) (m\a+ 

n,m 

h Vgo 2 e ,{(n - m+1) ^ +4}t (n| e lk2 ' r \m) \n) (m| a\ 



n.m 



h.c. 



(6) 



where we have introduced the dipole operators = 
12) (II) °2 = 12) (t|j and the quantities goi = A t oi'^oi/2fi, 
.902 = M02 ' E02/2/1. 

The large detuning condition allows the adiabatic elim- 
ination of the level |0) Under this condition, Ra- 
man transitions dominate the system. We also assume 
the system is cool enough to reach the Lamb-Dicke limit 
(kj> ■ r < 1) so we can expand 



e ll V 



1 + ikj/ x Xo {b + tf) 



(7) 



where f = 1, 2, xq = (h/2muj x ) 1 ^ 2 and m is the mass of 
the ion. Then the effective Hamiltonian in the interaction 
picture can be written 



i?/„ t = ^Ax |i> (i| + ftA 2 ir> <ri 
+ %6t|D(M + %*6||)<t| 



(8) 



where A; = \goi\ /A and g = ig^g^SkxX® / A with 
Sk x = k 2x - k lx . If we write S+ = |T) (||, S_ = ||> (T|, 
remove the terms for energy shifts and set g = g* , the ef- 
fective Hamiltonian has the anti Jaynes-Cummings form 

(9) 



HL f = hg (tfS+ + bS- 



This effective Hamiltonian (|9|) is the origin of the char- 
acteristic quantum dynamics (Rabi oscillations, collapses 
and revivals) of the system. Decoherence is the decay of 
the off-diagonal elements which represent the character- 
istic quantum dynamics, so we can use this Hamiltonian 
to explore some of the sources of decoherence in this in- 
teraction picture in the next section. 

When working in the interaction picture, it is con- 
venient to introduce the dressed states for the effective 
Hamiltonian (O): 



|p(n,l)> = -^(U,n) + |T,n + l)) J 
|p(n,2)> = -^(|4,n>-|T,n + l>) 

IT,o), 



(10) 

(11) 

(12) 

which are the eigenstates of the effective Hamiltonian. 
We write the eigenvalue of ( |l0| ) as E+, of (11) as E 1 ", of 
( |l2| ) as E 0l so we have E± = ±hg^n + 1 and E = 0. We 
write the reduced density operator in the dressed state 
basis as 

p 1 ® =£E<?(*)IpK/3)>M«,«)I 



+ 5>2o (*)lt,0) (<p(n,a)\ 

n,ct 

+ £p ( $(i)kp(m I /3)><T,0| 

m,{3 

+ poo (t)|T,o) (T,o| 



(13) 



for the boson quantum numbers n, m = 0, 1, 2, ... and the 
spin quantum numbers a, (3 — 1,2, where p"™, PqO'Po™' 
Poo are matrix elements. Then the population of the 
lower atomic state Pi is 



Pl(t) = l 1 



poo (f) + 2^Re[p^(t)] 



(14) 



Note that only the elements that are off-diagonal in terms 
of the spin quantum number (p™ (i)) and one diagonal 
element (poo (t)) contribute to P| in the dressed state ba- 
sis. Basically, the characteristic quantum dynamics ob- 
servable in the population of the lower state are due to 
the dynamics of elements that are off-diagonal in terms 
of the spin quantum number. 
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III. DECOHERENCE WITHOUT ENERGY 
RELAXATION 



A. Imperfect dipole transition 



We next consider the system with the effective Hamil- 
tonian (^|) described in the previous section now sur- 
rounded by the environment, that is, as an open system. 
Noise from the environment causes decoherence P,|l0|. 
We treat this open system by coupling to a quantum 
reservoir, which consists of an infinite number of many 
mode bosons 

H r = hJ2^iBjB l , (15) 
i 

where uji is the Ith reservoir frequency, Bj and Bi are 
the creation and annihilation operators of the reservoir 
bosons. Since the reservoir has infinitely greater degrees 
of freedom, the reservoir bosons are not affected by the 
system. Then the time evolution of the reservoir boson 
operators are given by 

B] (t) = e iHrt/H B J e -iH r t/H = e Mt B ^ ( lg ) 
Bi (t) = e iHrt/h Bie -iH r t/H = e -i"it Bh 

The system-reservoir coupling Hamiltonian is 

H sr = hJ2 C *J2 9si (b\ + B,) (18) 

where g s i is the coupling between a system operator C s 
and the Zth reservoir mode. The sum of the system op- 
erators C s has to be Hermitian. In the master equa- 
tion derived from the system-reservoir coupling ([l8|), the 
damping term consists of the system operators coupling 
to the reservoir operators. Thus the choice of the cou- 
pling between system operators and the reservoir deter- 
mines the effect of the reservoir. If we choose a system 
operator C s with the property 

C a \<p{n,a))=Y,Cfi\<P(n>P)) (19) 


the resulting master equation describes relaxation within 
the dressed states of the quantum number n, but not en- 
ergy relaxation between states with different n. This is 
because the time evolution of the density matrix elements 
in terms of \tp (n, a)) decouples for different n. The oper- 
ators S z , tfb are obviously of this type, as these operators 
do not even change the motional states \n) as well as the 
dressed state label n. The operator tfS + + bS- changes 
the motional state, but this operator does not change the 
dressed state occupation label n, so tfS+ + bS- is of this 
type, too. On the other hand, if we choose C s with 

C s \^(n,a))=Y l ^\f(m^n,P)), (20) 

then the resulting master equation describes transitions 
between states with different boson quantum numbers, 
which cause energy relaxation; S+ + 5_ and b + b^ are of 
this type. 



First, we treat the case when the system operator 
which is coupling to the reservoir is b'S+ + bS-. This 
case looks strange at first sight, but we can consider this 
as the result of imperfect dipole transitions between the 
level |0) and the level \j) (j —I, j) due to fluctuations of 
the driving laser intensity. We have previously described 
how phase fluctuations lead to decoherence and the de- 
struction of quantum revivals in the Jaynes-Cummings 
model [jn]]. This is one particular realization of "intrin- 
sic decoherence" in which off-diagonal density matrix el- 
ements relax without energy relaxation jll],[l2| ■ We note 
that these earlier results of ours apply to the experiments 
of Ref. § if the source of decoherence is relative phase 
fluctuations driving the ionic Raman transition. Here we 
analyse more general sources of decoherence. The im- 
perfect dipole transitions \\) 4=> |0) and |f) |0) are 
represented by 

of -^o-f+o-f 91 (Bj + Bt), (21) 
i 

4 - of + of ]T .9, (Bj + Bi) . (22) 
i 

We assume the system-reservoir coupling is weak enough, 
so we can neglect the terms that are second order in gi. 
Then we have, for example, 

S+ = -> a^af + a^a+ ^ 9l (bJ + B t ) . (23) 

I 

Thus the Hamiltonian describing the system-reservoir 
coupling is given by 

H sr = h {tfS+ + ^ g[ (B\ + B t ) (24) 
i 

where g[ = ggi . This Hamiltonian ( p4| ) can be interpreted 
as if the Rabi frequency (g) of the Jaynes-Cummings type 
system fluctuates due to the system-reservoir coupling as 

g^9 + J29l'( B l+ B l)- ( 25 ) 
/ 

Using a time convolution-less (TCL) formulation (this 
approach is described in detail in [^3| ) of the quantum 
damping theory || and the rotating wave approximation 
on the master equation Jl4[ the master equation in the 
interaction picture is 

§ i p I (t) = ±[H I eff ,p I (t)]+rp I (t) (26) 

where the damping term Bp 1 is given by |15j 
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rp 1 (t) = £ gf jT df I (<S/ (0 + (Bi (f ) B+) B ) 

x ( (-*') / (t) , + (-*') p 1 (t) , 6+5+] ) 

+ ((B i t (-t')5 i ) B + (B 1 (-i')5 | t > 5 ) 

x ( & t 5+,p I (t)bS-(-t / ) + bS-,p ! (t)tfS+(-t') ) 



(27) 



with 



tfS+ (t) = e lH 'ff t/h tfS+e-' lH 'tf t/R , 
bS-(t) = e* H 'ff t/h bS-e- zH 'ff f/h . 



(28) 
(29) 



The master equation (|26|) can be solved by expanding 
all system operators in terms of the dressed states under 
certain reservoir conditions. We require the reservoir to 
be the canonical distribution at temperature T and the 
time scale of the reservoir variables to be much shorter 
than the system variables so we can take the Markovian 
limit. 

We take the continuum limit of the reservoir modes, 



E 



dujD(uj) 



(30) 



and the damping function k (n) is 

K(n) ~ D(Ct n ) ■ g(£l n ), 



(35) 



which represents the effective contribution of the reser- 
voir bosons having frequency f2„. So the combination of 
these (k (n) {n (n) + 1/2}) represents the effective mean 
number of the reservoir bosons with frequency f2„. The 
quantity Ko^o is the contribution from zero frequency 
reservoir bosons. 

The analytical solution of the equation (p2|) is 



PIT (*) 



where 



-A„t 



{cos (B n t) - in n /B n sin (B n t)} p^ (0) 



fr J ^ t A n /B n sin (B n t)pft(p) 



(36) 



where D{u>) is the density of states of the reservoir. The 
corresponding continuum expression for gi is g(u>). The 
master equation is cast into a group of differential equa- 
tions for the density matrix elements. The time evolution 
of density matrix elements having different boson quan- 
tum numbers are decoupled due to the character of the 
coupling between the system operator and the reservoir. 
The time evolution of the diagonal elements (p^S) and 
the off-diagonal elements (p^, a ^ (3) having the same 
boson number (n) are also decoupled. 

To calculate the time evolution of P±, we only need 
the elements p™£ and poo- The equations for the time 
evolution of these elements are 



0_ 

m 
d 
m 



Poo (t) = 0, 



(31) 



Pl2 (*) 



-m nP ^ (t) 

(n + l){n(n) + 1/2} «(n)p^(t) 

2 (n + 1) K Q n p™ (t) 

(n + 1) {n(n) + 1/2} K (n)p^ (t) 



where 



(32) 



(33) 



n n = E+ - E~ = 2. 9 Vn+T. 
In (|3^) , the function of the reservoir bosons n (n) is 

n(n)=(e m ^ kBT -iy\ (34) 



A n = (n + 1) K (n) {n (n) + 1/2} + 2 (n + 1) n n 

= A* ipole , (37) 
B n = y/Q* - Aj. (38) 

With the chosen initial condition: | })(} | for the atom 
an d Pn\n)(n\ for the motional state, the real part of 
Eq. (WJ) is found to be 



Re[p^(t)} = 




jp- ) cos ( B.,l 



(39) 



where 9 n is a phase shift defined by 9 n = arctan (A n /B n ). 
Thus we see that the damping rate A n depends on the 
effective mean number of the reservoir bosons with fre- 
quency f2„, with a factor n + 1. The coupling to the 
reservoir also shifts the oscillation from H n t to B n t + 9 n . 
Since we assumed that the system-reservoir coupling is 
weak in our formulation, k (n) in A n must be much 
smaller than the Rabi frequency g. Thus we have relation 
A n <C B n ~ fl n . Under this condition, the population of 
the lower atomic state P{ is approximated to be 



n (*) = \ | 1 + E p» cos e ~ Ant } 



(40) 



which is in the same form as that seen in the experiments 

i. 
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B. Fluctuation of vibrational potential 

Next, we consider the case that the system couples to 
the reservoir via the system operator b'b. The system- 
reservoir coupling Hamiltonian is 



H sr = htfbj^g'i ( B \ + B > 



(41) 



This coupling describes fluctuations of the trap potential. 
Then the damping term is 

rp 1 (t) = ]T gf J* dt'{ ((Bj (t f ) Bi) b + (Bi (0 B}) B ) 

x [tfb (-if) p 1 (t) ,tfb] 

+ ({B}{-t')B l ) B + {B l (-t')B}) B ) 

x [tfb^ 1 (t)tfb(-t'j\ 1 (42) 



where 



b^b(t) = e lH 'ff t/r Vbe~ m 'ff t ^ 



(43) 



After expanding ( }42| ) in terms of the dressed states, 
the time evolution of and poo are 



—poo [t) = 0, 



(44) 



-- {n(n) + 1/2} K (n)p^(t) 

+ ± {n(n) + 1/2} n(n)p™(t). (45) 



The analytical solution of Re [p™] is given by (B9J) with 



A n = ~/s(n){n(n) + 1/2} = ^ 



(46) 



_B„ is defined by the equation (pq). We note that equa- 
tions (^J)-(^) coincide with those for the case of coupling 

to the reservoir via S z (H sr = HS Z ^ g\ (Bj + Bn). 



IV. ESTIMATION OF RESERVOIR VARIABLES 

The formulation of the decoherence rates A n lpole J37| ) 
and A n tb ( [46|) in the previous section shows that deco- 
herence originates in the relaxation of density matrix el- 
ements that are diagonal in the boson quantum num- 
ber but off-diagonal in the spin quantum numbers in the 
dressed state basis. The relaxation of the element p^fl 
for a ^ (3 is caused by the coupling to reservoir bosons 
at frequency of fi„ (= 2gy / n + 1). The effective contri- 
bution of reservoir bosons at frequency of fl n is therefore 
the key to understand the decoherence rate. 



The Rabi frequency g in the Boulder experiment 
is around 100 kHz, so reservoir bosons of order 100 kHz 
seem to be responsible for decoherence. These reservoir 
bosons have a much lower frequency than those responsi- 
ble for the case of spontaneous emission between internal 
atomic states, which is of order GHz here, and also popu- 
lation decay of motional states, which is of order 10 MHz. 
This low frequency nature of the reservoir bosons impor- 
tant here suggests that the reservoir may be at non-zero 
temperature whereas of course in the optical frequency 
regime, the reservoir is often approximated to be at a 
zero temperature. 

What are these reservoir bosons in the experiment? 
To consider this question and discuss the origin of deco- 
herence, we investigate the other characteristics, which 
the reservoir bosons should satisfy, by comparing our 
theoretical results to the experiment. For this purpose, 
we introduce normalised values, Af/ pole = A n lpole /g, 
A^ h = Afjg, u x = u x /g, f = k B T/hg, n„ = Q n /g, 
7o = 7o/<7j k (n) = k (n) / g. The normalised decoherence 
rates are 



A% p ° le = (n+l)k(n)f(n,f 
K ib = ~S(n) / (n,f) 



+ 2(n + l)K Q n , (47) 
(48) 



where 



1 



f[n,T)=n{n) + ^ = \ coth (49) 

The experimentally observed decoherence rate can be 
written as A n x = 70 (n + l) u . In the experiment, g/2n = 
94 kHz, 70 = 11.9 kHz, so we have 70 = 0.127/2tt. 

Let us further assume that k (n) described by ( |35| ) is 
given by a power d of the frequency Cl n Jt| , 



k (n) = afi„ = a (2\Ai + l) 



(50) 



where a is a damping constant (a <C 1). Some high fre- 
quency cut-off of the damping function is assumed to 
prevent divergence. These are the usual arguments given 
for the reservoir density of states. Generally, the case 
of d = 1 is known as the Ohmic case, since the choice 
of d gives a velocity dependent dissipation rate for the 
dissipative two-state system |fT|, and d = 3 is required to 
describe 3-D radiation fields | lq| . However, we do not 
restrict ourselves to d as an integer. We ignore the effect 
of zero frequency reservoir bosons: Ko«o = 0. 

The decoherence rates at n = have to coincide with 
70. So we have 70 = k(0)/(0,T) for the imperfect dipole 
transition case and 70 = k(0)/(0,T)/2 for the case of 
fluctuations of the vibrational potential. These condi- 
tions for 70 determine a when the value of 70 is given. 
Thus the decoherence rate is rewritten as 

Af pole = 70 (n + l) 1+d/2 / (n, f) // (0, f ) , (51) 

ijf = 70 (n + if 2 / (n, f) If (0, t) . (52) 
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The remaining unrestricted fitting parameters in our 
formulation are the normalised temperature T and 
the power dependency d in Eq.([34|). The value 

/ (o, fj // (n, T\ lies in the range 



(l + n)- 1/2 </(n,f) // (o,T) 



< 1. 



(53) 



We take the high temperature limit (T — > oo). 
This limit represents classical noise where the reser- 



( 



voir operators commute 
/ (n, Tj // (o, T\ becomes (n + 1 



BuBl 



= 

-1/2 



The value 



when T 



CO, 



so we have 



If^=7o(n + l) (d+1)/2 
if =7o(n + l) (d - 1)/2 . 



(54) 
(55) 



The linear form 70 (n + 1) is reached for d — 1 (the 
Ohmic case) for the imperfect dipole transition and for 
d = 3 (3-D radiation field) for the case of fluctuations of 
the vibrational potential. To get a power exponent of 0.7 
for n + 1 we need d ~ 0.4 for the imperfect dipole transi- 
tion case (Fig. ||) and d ~ 2.4 for the case of fluctuations 
of the vibrational potential. 



(1) (2) 




gt/2% qt/2% 

FIG. 2. The population of the lower atomic state P[ (t) against the normalised time gt/2n when (1) the initial internal state 
is ||) and the initial motional state is a Fock state with condition |1); (2) the initial motional state is a coherent state with 
condition \a = 3.0). For both figures, the dashed lines are for the case of no decoherence and the solid lines are for the case of 
an imperfect dipole transition with the coefficients d = 0.4 and 70 = 0.127/2-7T in the high temperature limit T — * 00. 



V. SUMMARY AND OPEN QUESTIONS 

In summary, we have shown here a model describing 
decoherence which destroys the characteristic quantum 
dynamics (collapse and revival) of the Jaynes-Cummings 
system without energy relaxation for the ion trap ex- 
periment The sources of decoherence are first intro- 
duced phenomenologically and then described by a mas- 
ter equation using a microscopic Hamilton formulation. 
We apply the model to the two possible actual sources of 
decoherence; one is the imperfect dipole transition, and 
the other is the fluctuation of vibrational potential. We 
solve the master equation under the Markovian approxi- 
mation and continuum limit of the reservoir modes. 

The analytical solution shows that decoherence is de- 
scribed by the reservoir bosons with frequency f2„ = 
Ig^fn + 1. Therefore, the effective contribution of the 
bosons at frequency of Q„ (which is of order 100 kHz) 
was found to be the key to understand the decoherence. 
This low frequency nature of the reservoir bosons com- 
pared to the spontaneous emission transition frequencies 
(which arc of order GHz) and population decay transi- 



tion frequencies (which is of order 10 MHz) suggests that 
the reservoir may be regarded to be at non-zero tempera- 
ture. If we assume the high temperature limit and certain 
density of states of the reservoir bosons, the decay rate 
coincides with that seen in the experiment 

To proceed further and investigate the origin of deco- 
herence for the Boulder ion trap experiment 0] , we would 
have to know a number of parameters: 

1 . The intensity fluctuations of the dye laser used for 
the stimulated Raman transition which seems to be 
order of 10 5 -10 6 Hz, so this may well be a candi- 
date for the fluctuation affecting the bosons. But 
we would need to know more about the frequency 
dependence of the intensity fluctuation around 100 
kHz to take the analysis much further. 

2. Noise from the trap potential from the radio fre- 
quency (100-200 kHz) radiation field is possible, 
but again we would need to estimate the density 
of states for this case to be more precise. 



G 



3. Possibility of quantum noise (noise at finite T) re- 
mains a potential candidate to explain this deco- 
herence. 

We defer further consideration of all these until the un- 
derlying parameters are better understood. Since this 
paper was submitted for publication we learnt of related 
work by Schneider and Milburn |17| , and by James mM . 
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